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1. INTRODUCTION. In 0, Euler analyzed the Riccati equation to prove that the 
number e has the continued fraction expansion 



e= [2,1,2,1,1,4,1,1,6,1,1,8,...] =2 



1 



The pattern becomes more elegant if one replaces the initial 2 with 1,0, 1 , which yields 
the equivalent continued fraction 

[1,0,1,1,2,1,1,4,1,1,6,1,1,8,1,1,10,1,...] (1) 

because 

1 + —. =2 + .-.. 



One of the most interesting proofs is due to Hermite; it arose as a byproduct of his 
proof of the transcendence of e in J5|. (See |6| for an exposition by Olds.) The purpose 
of this note is to present an especially short and direct variant of Hermite's proof and 
to explain some of the motivation behind it. 

Consider any continued fraction [ao,ai,02, ■ ■ ■]■ Its ith convergent is defined to 
be the continued fraction [ao, a\, . . . , a,]. One of the most fundamental facts about 
continued fractions is that the ith convergent equals p-,/q„ where p, and q- t can be 
calculated recursively using 

Pn = GnPn-1 + Pn-2, In = d„q n -\ + q n -2, 

starting from the initial conditions po = ao, #o — 1> Pi = ^ofli + 1, and q\ — a\. For 
a proof see |4,, p. 130] or any introduction to continued fractions. 
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2. PROOF OF THE EXPANSION. Let [a , a u a 2 , . . .] be the continued frac- 
tion Q. In other words, 03,+! = 2i and = an+% = 1, so p, and g,- are as follows: 



i 





1 


2 


3 


4 


5 


6 


7 


8 


Pi 


1 


1 


2 


3 


8 


11 


19 


87 


106 


<li 


1 





1 


1 


3 


4 


7 


32 


39 



(Note that q x = so />i/gi is undefined, but that will not be a problem.) Then />, and 
q, satisfy the recurrence relations 

Pin = Pin-\ + Pin-2, <?3n = #3n-l + <?3n-2) 

P3«+l = 2nj? 3 „ + p 3 «-l, <?3»+l = 2«^ 3 « + <?3;,-l, 

P3m+2 = /?3n+l + Pin, lin+2 = <?3n+l + l3n- 

To verify that the continued fraction Q equals e, we must prove that 

lim — = e. 

Define the integrals 

.1 v"C 

A„ 

B„ = I ^— —e x dx, 

Jo n\ 

c n =r^ x -p n+1 e * d , 





—e 


n\ 




x" +1 (x- 


1)" 


n\ 




x n (x-l 


!«+! 



Proposition 1. For n > 0, A„ = q 3 „e — p^„, B n = p^ n+ \ — q3„+ie, and C n = P3„+2 — 
qin+ie- 

Proof. In light of the recurrence relations cited earlier, we need only verify the initial 
conditions Aq = e — 1 , Bq = 1 , and Co = 2 — e (which are easy to check) and prove 
the three recurrence relations 

A„ = —B„-i — C„_i, (2) 
B n = -2nA„+C„- U (3) 
C n =B„-A„. (4) 

Of course, © is trivial. To prove (jjjl (i.e., A„ + 5„_i + C„_i = 0) integrate both 
sides of 

x"{x-l)" x x"(x-l)"-' x x"- l (x-l)" x _ d (x"(x-l)" x 
n\ 6 + (n-1)! 6 H («-l)! & ~ dx \ n\ & 

which follows immediately from the product rule for derivatives. To prove Q (i.e., 
B„ + 2nA n — C„_i = 0) integrate both sides of 

x" +l (x-l)" , x n (x-\) n x n -\x-lf r d /x"(x-l)" +l 

-e x + In— -e x — -^-e x = — 1 K 

n\ n\ in — 1 ) ! dx 
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which follows from the product rule and some additional manipulation. This completes 
the proof. ■ 



The recurrences (|2j and © can also be proved by integration by parts. 
Theorem 1. e = [1,0, 1, 1,2, 1, 1,4, 1, 1,6, 1, 1,8, 1, 1, 10, 1, . . .]. 
Proof. Clearly A n , B„, and C„ tend to as n — > °°. It follows from Proposition^that 



3. MOTIVATION. The most surprising aspect of this proof is the integral formulas, 
which have no apparent motivation. The difficulty is that the machinery that led to 
them has been removed from the final proof. Hermite wrote down the integrals while 
studying Pade approximants, a context in which it is easier to see how one might think 
of them. 

Pade approximants are certain rational function approximations to a power series. 
They are named after Hermite's student Pade, in whose 1892 thesis [7] they were 
studied systematically. However, for the special case of the exponential function they 
are implicit in Hermite's 1873 paper |5| on the transcendence of e. 

We will focus on the power series 




since it is the relevant one for our purposes. A Pade approximant to e : of type (m, n) is 
a rational function p(z)/q(z) with p(z) and q{z) polynomials such that deg p{z) < m, 
degq(z) < n, and 



as z — > 0. In other words, the first m + n + 1 coefficients in the Taylor series of 
p(z)/q{z) agree with those for e z . One cannot expect more agreement than that: p(z) 
has m + 1 coefficients and q(z) has n + 1, so there are m + n + 2 degrees of free- 
dom in toto, one of which is lost because p{z) and q{z) can be scaled by the same 
factor without changing their ratio. Thus, one expects to be able to match m + n + 1 
coefficients. (Of course this argument is not rigorous.) 

It is easy to see that there can be only one Pade approximant of type (m,n): if 
r(z) /s(z) is another, then 



lim qie — p, 



= 0. 



Because 



q, > 1 when i > 2, we see that 



e = lim^ = [1,0,1,1,2,1,1,4,1,...] 

i->°° (7; 



as desired. 




p{z)s(z) - q{z)r{z) _ p(z) 



q(z)s(z) q(z) 




0(z m+n+l ) . 
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Because p(z)s(z) — q(z)r{z) vanishes to order m + n + 1 at z = but has degree 
at most m + n (assuming degr(z) < m and degs(z) < n), it must vanish identically. 
Thus, p(z)/q(z) = r(z)/s(z). 

We have no need to deal with the existence of Pade approximants here, because it 
will follow from a later argument. One can compute the Pade approximants of e z by 
solving simultaneous linear equations to determine the coefficients of the numerator 
and denominator, after normalizing so the constant term of the denominator is 1 . 

The usefulness of Pade approximants lies in the fact that they provide a powerful 
way to approximate a power series. Those of type (m, 0) are simply the partial sums 
of the series, and the others are equally natural approximations. If one is interested in 
approximating e, it makes sense to plug z = 1 into the Pade approximants for e z and 
see what happens. 

Let r m , n {z) denote the Pade approximant of type (m, n) for e z . Computing continued 
fractions reveals that 



'1,1 (1) = 


[2,1], 


'1,2(1) = 


[2,1,2], 


'2.1 (i) = 


[2,1,2,1], 


'2,2(1) = 


[2,1,2,1,1], 


'2,3(1) = 


[2,1,2,1,1,4], 


'3,2(1) = 


[2,1,2,1,1,4,1], 


'3,3(1) = 


[2,1,2,1,1,4,1,1 



etc. In other words, when we set z = 1, the Pade approximants of types (n,n), (n, 
n + 1), and (n + l,n) appear to give the convergents to the continued fraction of e. 
There is no reason to think Hermite approached the problem quite this way, but his 
paper |5| does include some numerical calculations of approximations to e, and it is 
plausible that his strategy was informed by patterns he observed in the numbers. 

It is not clear how to prove this numerical pattern directly from the definitions. How- 
ever, Hermite found an ingenious way to derive the Pade approximants from integrals, 
which can be used to prove it. In fact, it will follow easily from Proposition^ 

It is helpful to reformulate the definition as follows. For a Pade approximant of 
type (m,n), we are looking for polynomials p{z) and q{z) of degrees at most m and 
n, respectively, such that q{z)e z — p(z) = O (z m+ " +1 ) as z — > 0. In other words, the 
function 

q(z)e z - p(z) 
z >— > 

-,m+«+ 1 

must be holomorphic. (Here that is equivalent to being bounded for z near 0. No com- 
plex analysis is needed in this article.) 

One way to recognize a function as being holomorphic is to write it as a suitable 
integral. For example, because 

— — + 1 = / xe ZJC dx, 
z l Jo 

it is clear that z >— > ((z — l)e z + l)/z 2 is holomorphic. Of course, that is unnecessary 
for such a simple function, but Hermite realized that this technique was quite powerful. 
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It is not clear how he thought of it, but everyone who knows calculus has integrated an 
exponential times a polynomial, and one can imagine he simply remembered that the 
answer has exactly the form we seek. 

Lemma 1. Let r(x) be a polynomial of degree k. Then there are polynomials q{z) and 
p(z) of degree at most k such that 

1 i \ zx j l(z)e z ~ P(z) 
r(x)e dx = — 



o ' ' z k+l 



Specifically, 



and 



q ( z ) = r (l)z* - r\l)z k - x + r"(l)z k - 2 



p(z) = r(0)z k - r'(0)z k - 1 + r"(0)z k - 2 
Proof. Integration by parts implies that 

/ r{x)e~ x dx= — — / r'(x)e <x dx, 

Jo z z Jo 

from which the desired result follows by induction. ■ 

To get a Pade approximant p{z)/q(z) of type (m,n), we want polynomials p(z) 
and q(z) of degrees m and n, respectively, such that 

q(z)e z - p(z) 



is holomorphic. That suggests we should take k = m + n in the lemma. However, if 
r(x) is not chosen carefully, then the degrees of q(z) and p(z) will be too high. To 
choose r(x), we examine the explicit formulas 

q( z ) = r{l)z m+n - r'(l)z"' + "- 1 + r"(l)z ffl+ "- 2 

and 

p( z ) = r(0)z" ,+ " - r'(0)z" ,+ "- 1 + r"(0)f 1+ "- 2 . 

The condition that degg(z) < n simply means r{x) has a root of order m at x = 1, and 
similarly deg p(z) < m means r{x) has a root of order n at x = 0. Since degr(x) = 
m + n, our only choice (up to a constant factor) is to take 

r(x) =x n {x-\)"\ 

and that polynomial works. Thus, 

q{z)e z - p{z) 



[ x"\ 
Jo 



x-l)"'e zx dx 



ym+n+ 1 ' 



where p(z)/q(z) is the Pade approximant of type (m, n) to e : 
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Setting z = 1 recovers the integrals used in the proof of the continued fraction 
expansion of e, except for the factor of l/n\, which simply makes the answer prettier. 
Fundamentally, the reason why the factorial appears is that 



has integral coefficients. Note also that up to a change of variables, this expression is 
the Rodrigues formula for the Legendre polynomial of degree n Q] p. 99]. 

Natural generalizations of these integrals play a fundamental role in Hermite's proof 
of the transcendence of e. See HJ p. 4] for an especially short version of the proof or 
chapter 20 of |8] for a more leisurely account (although the integrals used there are 
slightly different from those in this paper). 
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